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Abstract. We investigate spherically symmetric solutions with pressure and discuss the existence 
of a dividing shell separating expanding and collapsing regions. We perform a 3+1 splitting and 
obtain gauge invariant conditions relating not only the intrinsic spatial curvature of the shells to the 
ADM mass, but also a function of the pressure which we introduce that generalises the Tolman- 
Oppenheimer-Volkoff equilibrium condition. We consider the particular case of a Lemaitre-Tolman 
dust models with a cosmological constant (a A-CDM model) as an example of our results. 
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INTRODUCTION 

Cosmological formation of structure assumes the collapse of inhomogeneities, via grav- 
itational instability, into "bound" structures, with the underlying idea that they depart 
from the cosmological expansion. This approach usually models overdensities with 
closed patches embedded in Friedman backgrounds, in particular in the spherical col- 
lapse included in the Press & Schechter scheme [1]. Birkhoff's theorem is invoked to 
justify their independent evolutions [2], while it is proved for asymptotically flat space- 
times [3]. Such separation is reminiscent of the concept of trapped surfaces [4]. 

In this work we define general conditions for the existence of a shell separating 
contraction from expansion, and will illustrate our results with a simple example of 
perturbed A-CDM models. We adopt the Generalised Painleve-GuUstrand (hereafter 
GPG) formalism used in Lasky & Lun [5], which involves a 3 + 1 splitting (ADM) 
and the consideration of gauge invariants kinematic quantities [6] . 

ADM APPROACH TO LTB MODELS IN GPG SYSTEM 

We consider a spherically symmetric Generalised Lemaitre-Tolman-Bondi (Generalised 
LTB, or GLTB) metric to include pressure. Performing an ADM 3-1-1 splitting in the 
GPG coordinates [5] , the metric reads 

ds^ = -a{t,rfdt^ + \ — -m{t,r)dt + drf + ?d^^ , (1) 
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characterised where a is a lapse function, /3 is a shift, and £ is a curvature-energy 
characterising the curvature of the spatial surfaces orthogonal to the direction = 
(—oc, 0,0,0) of the flow. For a perfect fluid, the projected Bianchi identities T^.^ = 




yield the energy density conservation equation when we project along the flow n^, and 
the Euler equation when we project orthogonally to the latter. Using the projection h^, 
and denoting the density Lie derivative along the flow =Sf„p, we have: 

n%% = - if„p - (p +P) = 0, I^T^:^ =0 ^ i" = - (p +/') - , (2) 

where P is the pressure, the radial derivatives are denoted by a prime, the time 
derivatives are represented by a dot, and the expansion is 0. 
Introducing the ADM (also called Misner-Sharp) mass 

M = r^{l+E){\nay -47tPr^ + -Ar^ + r^^n(^j , (3) 

where A is the cosmological constant, we write the Einstein's field equations (EFEs) as 
Lie derivatives along the flow 



(4) 
(5) 
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The system becomes then closed when an equation of state (EoS) f{p,P) = is 
supplied. The A term can be absorbed and many fluids mass equations written with 

each component using the ^ term for the overall sum of the masses. 

SEPARATING COLLAPSE FROM EXPANSION 

To characterise the separation of expansion from collapse, we prefer to use here GLTB 
coordinates in which the flow direction reduces to time by choosing /3 = — r. We thus 
have 

ds^ = -a{T,Rf {drtfdT^ + + rW, (7) 

M=l54nPr^,Er' = 2I5^-^P', (--^ =E + 2— + ^Ar^. (8) 

p ~t~ P \^ CC J r 3 

There are two situations one should consider in parallel. On the one hand, we look for 

the vanishing of the gauge invariant expansion 0, defined as n"-a, since our goal is to 



separate an inner collapsing, spherical region from the outer expanding universe. On the 
other hand, the total ADM mass of the spherical region that departs from the expansion 
flow should be conserved. This is indeed suggested by the dust case where that happens 
for every shell. 

We denote with ★ the dividing shell. The analysis of the kinematic quantities reveals 
that is linked to the shear a, in that at the ★ shell 



= -3 
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(= when ^„M{t, r^{t)) = 0) . (9) 



(10) 



a + ^-|<^ 0,, + 3a,, = -3- 
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Moreover the generalized Friedman constraint 

(3)^ + ^02 =6fl2 + I6np + 2A 

tells us that the vanishing of only happens in region of positive 3-curvature ^^^R. On the 
other hand, if we demand that the separating shell has a dust-like vanishing mass/energy 
flow, i.e., has a conserved ADM mass along n'': ^ 

yt,^nM{t,r^{t))=0 ^yt,E = -2— <0. (11) 

We remark that M refers to the total ADM mass, thus including a cosmological constant. 

Taking into account that the equilibrium of static spherical configurations requires 
the satisfaction of the Tolman-Oppenheimer-Volkoff equation of state [7], we are led to 
define a generalized TOV function 

T 

a 



TOV 



1+E. M 1 

-F +47tPr+^--Ar 



(12) 



_p+P' ■ ■ 3 

which reduces to the usual TOV equation when it vanishes. 

The radial behaviour of the 7»c-shell is then similar to a turnaround shell. Indeed 
= — leads to a null radial velocity = while its acceleration reveals the 
importance of the TOV parameter 

TOV^l 

' ; ™,* = -a^TOV^. (13) 

The local staticity of the ★-shell is then shown to be equivalent to having a local TOV 
equation on this limit shell r . m " 
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AN EXAMPLE: ACDM 

A simple illustration of our result is given by the case of dust with a A. There is then 
no pressure gradients and M^i^st and E are conserved, which simplify the analysis (i.e. 
a — 1), and allow us to perform a kinematic analysis (see Fig. 1) per shell of Eq. (8-c). 
An effective potential is defined by V{r) = — ^ — ^r^, so that when r = we have 
indeed V{r) — E: 
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For each shell there is a virtual static state at ru, 
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only depends on A and M{R) . We also recall that, in this model, TOV 

We then need to choose initial conditions p,, which sets the En^ profile, and Vi, which 
sets the Ei profile. The intersection of En^ with Ej will be static at turnaround, hence 
defining our limit shell. We use two sets of cosmologically motivated initial conditions: 
one with p, as an NFW profile [8] and Ei as a parabola, and the second, with a cuspless 
power law as p, and a Hubble flow for v;. 



CONCLUSIONS 

Using non- singular, Generalized Painleve-GuUstrand coordinate formulation of the 
ADM spherically symmetric, perfect fluid system [5] we have shown [9] that the exis- 
tence of shells locally separating between inner collapsing and outer expanding regions, 
is governed by the condition that the combination of expansion scalar and shear 6 +3a 
should vanish on the shell. The ADM mass of the shell is then conserved. This condition 
requires that the separating shell must be located in an elliptic (E < 0) region. Moreover, 
for that shell to exist for certain over time, we have shown that the TOV equation must 
be locally satisfied. We argue that this local condition is global in a cosmological context 
(FLRW match at radial asymptote). Given appropriate initial conditions, this translates 
into global separations between an expanding outer region and an eventually collapsing 
inner region. We present simple but physically interesting illustrations of the results, a 
model of Lemaitre-Tolman dust with A representing spherical perturbations in a ACDM 
with two different initial sets of cosmologically interesting conditions consistent with 
known phenomenological constraints [9, and Refs. therein]: an NFW density profile 
with a simple curvature profile going from bound to unbound conditions, and a non 
cuspy power law fluctuation with initial Hubble flow. We show, for these models, the 
existence of a global separation. We argue that these shells are trapped matter surfaces 
[9] and that they constitute the validity locus to an analog to Birkhoff 's theorem. 
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